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Preface 


Classical  plate  analysis  generally  assumes  initially 
perfect  plates  when  considering  collapse  due  to  inplane 
stresses.  This  report  incorporates  the  assumption  of  initial 
imperfections  in  geometry  in  an  analytical  investigation 
of  collapse  involving  elastic,  elastic-plastic  and  ortho¬ 
tropic  plates. 

I  wish  to  express  my  sincere  appreciation  to  Dr. 
Anthony  Palazotto  for  his  valuable  instruction,  encouragement 
and  timely  advice.  1  particularly  want  to  thank  my  wife 
and  family  for  their  patience  and  understanding  during  my 
preoccupation  with  this  thesis. 
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Abstract 


An  analysis  is  made  of  the  effect  of  initial  imper¬ 
fections  in  geometry  on  the  collapse  of  simply-supported 
rectangular  plates.  A  Reissner-type  variational  principle 
is  employed  to  evaluate  the  load  versus  lateral  displacement 
of  elastic,  elastic-plastic  and  orthotropic  plates  subject 
to  compressive  loading  along  two  edges.  Imperfection  size 
and  shape,  plate  thickness  and  aspect  ratio,  and  for  ortho¬ 
tropic  material,  the  ply  orientation  is  examined.  The 
results  indicate  that  the  load-displacement  curve  is 
sensitive  to  the  initial  imperfection,  particularly  for 
elastic-plastic  materials.  The  effect  of  plate  size  and 
aspect  ratio  are  similar  for  all  the  plates  considered 
with  the  collapse  stresses  higher  for  elastic  plates  than 
for  elastic-plastic  plates.  The  effect  of  plasticity 
becomes  insignificant  for  thin  plates  when  compared  with 
plates  with  elastic  properties.  In  general,  it  can  be 
stated  that  the  imperfection  function  determines  the  plate's 
collapse  load.  For  orthotropic  plates  the  effect  of  ply 
orientation  is  significant  at  aspect  ratios  less  than  1 
but  relatively  insignificant  at  higher  values  when  comparing 
the  minimum  collapse  stresses.  The  difference  is  dependent 
upon  the  magnitude  of  the  ratio  of  Young's  moduli  in  the 
major  material  axes. 
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EFFECT  OF  IMPERFECTIONS  ON  THE  COLLAPSE 


OF  RECTANGULAR  PLATES 
USING  VARIATIONAL  CALCULUS 


I,  Introduction 


Background 

The  first  problems  in  elastic  instabi  ty  were  solved 
over  200  years  ago  by  Euler  (1)  working  with  the  lateral 
buckling  of  compressed  members.  Since  that  time  the  science 
of  structural  stability  has  received  attention  by  numerous 
individuals.  Because  of  its  application  as  a  structural 
member  in  modern  structures,  the  rectangular  plate  is  of 
particular  interest  in  this  thesis. 

In  the  late  1800's  Bryan  (2)  studied  plates,  applying 
the  theory  to  the  buckling  of  the  sides  of  ships.  In  the 
early  1800's  Timoshenko  (3)  treated  numerous  types  of 
loading  and  boundary  conditions  utilizing  both  the  energy 
approach  and  the  solution  of  the  differential  equations. 

In  his  1940  work.  Hill  (4)  constructed  charts  describing  the 
critical  compressive  stresses  for  flat  rectangular  elastic 
plates.  The  results  of  the  above  cited  work  compare  favor¬ 
ably  with  experimentation. 

In  1942  Lundquist  and  Stowell  (5)  presented  the  first 
unified  treatment  of  the  elastic-plastic  compressive  buckling 
problem,  using  both  the  differential  equations  and  energy 
methods.  However,  much  of  the  classical  work  with  elastic- 
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plastic  plates  resulted  in  critical  compressive  stresses  ! 

which  were  significantly  higher  than  experimental  results. 

It  was  shown  ioy  Onat  and  Drucker  (6  7)  for  the  axially 

compressed  cruciform  column  and  by  Neale  (8)  for  the  ! 

i 

i 

cylindrical  shell  under  torsion  that  these  theoretical  j 

t 

stresses  could  be  brought  much  closer  to  experimental  j 

results  by  including  in  the  analysis  initial  geometric  'j 

imperfections.  Consequently,  the  idealization  of  a  geo¬ 
metrically  perfect  system  adopted  in  the  classical  approach 
no  longer  remains  acceptable.  In  his  work  presented  in  i 

1974,  Neale  (9)  employed  the  variational  principle  and  a 
Ritz  method  incorporating  imperfections  into  the  analysis 
of  the  plastic  buckling  of  rectangular  plates.  His  results 
show  much  closer  agreement  to  experimental  data. 

The  1960 's  and  1970 's  have  seen  the  increasing  use  of 

composite  material  particula*.ly  in  structures  requiring  high  I 

1 

strength-to-weight  ratios.  These  composite  materials  often  I 

i  ; 

possess  orthotropic  properties,  and  consequently  modified 
theories  for  the  solution  of  problems  where  composite  mate-  i 

rials  are  used,  are  required.  Ashton  and  Whitney  (10)  present 
the  theoretical  development  and  solutions  for  plates  fabri¬ 
cated  of  thin  layers  of  anisotropic  material.  Jones  (11)  , 

and  Chamis  (12)  present  the  major  aspects  of  the  science  and 
technology  of  composite  materials  that  are  the  basics  of 
todays  technology. 

These  three  types  of  material,  elastic,  elastic-plastic 
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and  orthotropic,  form  the  basis  of  the  spectrum  of  material 
properties.  With  the  exception  of  Neale's  work  with  elastic- 
plastic  plates,  little  has  been  published  that  illustrates 
the  effect  of  initial  geometric  imperfections  on  the  collapse 
of  rectangular  plates.  Consequently,  it  is  this  subject  to 
which  this  thesis  is  devoted. 

Purpose 

The  purpose  of  this  thesis  is  to  examine  analytically, 
the  effect  of  small  and  consequently  unavoidable  imperfections 
in  geometry  on  the  collapse  of  simply  supported  rectangular 
plates  subject  to  in-plane  compression  along  two  edges. 

Three  different  types  of  material,  elastic,  elastic-plastic, 
and  orthotropic,  will  be  compared  and  the  effect  of  aspect 
ratio  (i.e.  length/width),  length/ thickness ,  shape  of  initial 
imperfections,  size  of  imperfection  and  for  the  orthotropic 
case,  ply  orientation,  will  be  examined. 

The  analytical  technique  that  will  be  used  consists 
of  the  formulation  of  a  "rate"  or  incremental  boundary-value 
problem,  which  can  be  solved  in  an  approximate  manner  using 
a  variational  principle  (13)  in  conjunction  with  the  method 
of  Ritz.  A  stepwise  solution  of  this  rate  problem  furnishes 
the  pertinent  load-deformation  behavior  (9) .  Using  this 
formulation  the  problem  reduces  to  a  set  of  five  algebraic 
equations  which  can  be  solved  incrementally  using  a  minimum 
of  computer  time  or  possibly  an  advanced  hand  calculator  or 
desk  computer. 
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Scope 

All  the  plates  examined  are  0.1  inches  thick.  Length- 
to-thickness  ratios  vary  from  10  to  400  and  aspect  ratios 
vary  from  0.5  to  4.0.  For  each  case  an  initial  imperfection 
of  lO”®,  10”**,  10”^  inches  is  evaluated.  The  effect  of  the 
shape  of  the  imperfection  is  evaluated  using  m  and  n  of 
1,  3,  5,  where  m  is  the  number  of  half  sine  waves  in  the 
x-direction  and  n  is  the  number  of  half  sine  waves  in  the 
y-Kiirection.  For  further  explanation  see  results  and 
discussion.  For  the  orthotropic  case,  a  single  ply  with 
fiber  orientation  at  0°  and  90°  with  respect  to  the  applied 
load  is  evaluated. 

Computer  values  of  critical  load  are  compared  with 
Timoshenko  (3)  for  the  elastic  case,  Neale  (9)  for  the 
elastic-plastic  case  and  Chamis  (14)  ,  and  Mandell  (15)  for 
the  orthotropic  case. 
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II .  Theory 

Since  the  case  of  elastic-plastic  material  is  the 
most  general  of  the  three  materials ,  the  theoretical  anal¬ 
ysis  will  center  on  plasticity  with  appropriate  modification 
made  for  the  elastic  and  orthotropic  cases. 

Reissner  (16  -  19)  provided  the  starting  point  for 
the  derivation  with  his  general  variational  theorem  for 
stresses  and  displacements.  Reissner 's  theorem  presents 
the  equilibrium  equations  and  stress-strain  relations  as 
Euler  Equations  of  the  variational  problem.  Boundary 
conditions  for  stresses  and  displacements  are  the  natural 
boundary  conditions .  The  method  permits  independent 
selection  of  both  the  stress  distributions  and  displacement 
quantities  in  the  approximate  solution  to  the  boundary- 
value  problem. 

Since  the  undeformed  configuration  of  the  material 
is  often  used  as  the  reference  frame,  the  Lagrangian  form 
of  the  equations  is  more  appropriate  than  the  Eulerian  form, 
which  refers  to  the  deformed  configuration.  In  conjunction 
with  this,  suitable  definitions  of  stress  and  strain  are  the 
Kirchhoff's  stress  tensor  and  Green's  strain  tensor  (20). 

General  expressions  for  the  Euler  Equations  of  the 
variational  principal  and  natural  boundary  conditions  can 
be  developed  using  the  aforementioned  relations.  A  complete 
derivation  can  be  found  in  reference  (13) .  These  relations 
are  subsequently  applied  to  rectangular  plates.  In  accord- 
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ance  with  Neale's  work  (9),  rectangular  cartesian  coordinates 
X,  y  describe  the  location  of  points  on  the  midplane  of  the 
plate  and  z  describes  the  distance  from  the  midplane,  as 
shown  in  figure  1.  By  assuming  normal  plane  surfaces  remain 
plane,  the  displacement  components  are  specified  as  follows: 

Ujj  =  u  -  23w 
dx 

u  =  V  -  z3w  (1) 

y  3y 

Ug  =  w 

where  u,  v,  w  are  the  displacement  components  in  the  x,  y,  z 
directions  respectively.  The  Green  strains  associated  with 
these  displacements  are  (20) : 

Ejjx  =  ^  -  z3^vsf  + 

3x  ¥x^'  \3x/ 

=  3v  -  z  3^w  + 

3y  Kdy) 

Ejjy  =  h  ^  ^  -  2z  3  ^w  +  ^  ^1  (2) 

3y  3x  3x3y  3x  3y  j 

The  products  and  squares  of  the  derivatives  incorporating 
the  transverse  displacement  w,  are  the  only  nonlinear  terms 
retained.  The  strain  rates  can  be  determined  from  equation 
(2)  yielding  results  which  assume  homogeneity  in  time. 


XX 


3u  _  3^w  3w  3w 
3x  "  ^ax®  "Sx  3x 


Eyy 


_  3^w  .  3w  ^ 

9y  3y^  3y 


xy 


L  f3_u  ,  3v  27  3w  3w  3w  3w 

^  'Sy  ^  ~  3x3y  3x  3y  3x  3y 


(3) 


where  the  dot  above  the  term,  (*),  denotes  an  incremental 
function.  The  incremental  stress  and  strain  constitutive 
equations  may  be  stated  by  the  following  matrix  expressions; 

*“11  *“i  a  *“13 

Cj2  C22  C2 3 

*“13  ^2  3  ^33 
2  2 

in  which  the  quantities  t  t  are  considered 

xz  f  y z  ^  z  z 

negligible  by  the  assumption  of  plane  stress. 

For  an  isotropic  solid  characterized  by  elastic- 
plastic  response  incorporating  the  Von  Mises  incremental 
theory,  one  may  obtain  the  •  coefficients  in  Eq. (4)  (21) : 
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Cs  3 


1+V 

E 


+  2Gt 


(5) 


where  E  is  Young's  modulus  and  v  is  Poisson's  ratio.  If  it 
is  assumed  that  the  uniaxial  stress-strain  can  be  described 
by  a  Ramberg-Osgood  relation  of  the  form 


E 

XX 


Txx 


(6) 


where  E^  and  k  are  material  constants,  then  for  plastic 
loading,  the  hardening  function  G  can  be  written  as 


G 


(7) 


The  quantity  Jj  is  stated  as  (21) , 


J,  =  1/3(t^^^-t 


XX 


XX  yy 


+  T 


yy 


)  +T 


xy 


(8) 


To  simplify  the  analysis,  k  is  taken  to  be  3  so  that  the 

hardening  function  becomes 

G  =  27/4E  ^  (9) 

0 


which  remains  constant  throughout  the  load  application. 

It  may  be  observed  that  the  strain  increments  in  Eq(4)  are 
functions  of  the  known  stress  components  and  a  set  of  new 
incremental  stresses.  This  set  of  constitutive  relationships 
is  known  as  the  Prandtl-Reuss  flow  law  of  plasticity. 

Likewise  when  evaluating  an  isotropic,  elastic  material 
under  plane  stress,  matrix  becomes  (22); 
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Cjj  ~  g  ,  C22  £  »  ^2  3  “  ^ 

„  1  +  V 

C3  3  =  “E"  (10) 

As  an  added  feature  to  the  use  of  Eq{4),  one  may 
modify  the  expressions  in  order  to  account  for  orthotropic 
material  properties  associated  with  composites.  Using  6 
as  the  angle  between  the  major  material  axis  and  the  x-axis, 
the  coefficients  may  be  stated  as  (11) : 


Cji  =  Sjj  cos‘*0  +  (2Sj2  +  Sgg)  sin^0  cos^0  +  S22  sin‘*0 

Cj2  =  Sj2(sin'*0  +  cos'*©)  +  (Sjj  +  S22  ~  Sgg)  sin^0  cos^0 

Ci3  =  (2Sn  -  2Si2  -  Sgg)  sin0  cos’0  -(2S^^  -  -  Sgg) 

^sin^0  cos©) 

C22  =  Sjj  sin**©  +  (2Sj2  +  ^66^  sin^©  cos^©  +  5^2  cos'*© 

C23  =(2S^j  -  25^2  -  Sgg)  sin^©  cos©  -  (2S22  -  25,2  "  S^g) 

( sin©  cos ^©  ) 

C33  =  (2Sjj  +  2S22  ~  ^^12  ~  sin  6  cos  ©  +  Sgg 

(sin"*©  +  cos'*©)  (11) 

where 
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( 


( 


( 


The  quantities  and  are  the  Young's  moduli  in  the 
longitudinal  and  transverse  material  directions.  The  term 
is  Poisson's  ratio  for  transverse  strain  in  the  j- 
direction  when  stressed  in  the  i-direction.  G, ,  is  the 
shear  modulus  in  the  1-2  plane. 

Thus,  Eq.  (4)  is  useful  for  any  general  material 
as  long  as  some  constitutive  relationship  can  be  found. 


Formulation  of  Rate  Boundary-Value  Problem  for  Plates 

Applying  the  general  solution  of  the  rate  boundary- 
value  problem,  the  functional  1°  can  be  stated  for  a  simply 
supported  rectangular  plate  as  (9) 
t/2  b  a 


1°  = 


-t/2  -b  -a 


^  ""ij  ^k,j 

(13) 


where  t  -  .  E.  ,  =  t  E  +  2t  E„„  +  t  E 

13  ij  XX  XX  xy  xy  yy  yy 

■^ii  '^k  i  ^k  i  "  "^xx  ^"^xv  +  lyy  ^  ^ 

1]  K,1  K,3  XX  xy  g^  YY  g^ 

2  W  (t)  =  Cjj  +  2Cj2  Tjjjj  lyy  ^ 2  2  '^yy  ^^13  ^XX  ^xy 


'2  3  xy  yy  ''3  3  xy 


(14) 


with  specified  in  equations  (5) ,  (6) ,  (11)  for  elastic- 
plastic,  elastic  and  orthotropic  materials  respectively. 

The  general  variational  principal  states  the  functional 
1°  is  stationary  in  the  vicinity  of  the  solution.  The  boundary 
value  problem  consists  of  determining  the  incremental  stress 


11 


«  «  * 


and  displacement  functions,  u,  v,  v,  which 

satisfy  the  boundary  conditions  corresponding  to  the  specified 
edge  displacements.  Using  the  method  of  Ritz,  an  approximate 
solution  for  these  incremental  quantities  can  be  obtained. 

The  problem  considered  is  one  in  which  a  u  displacement 
is  initiated  in  a  plate  characterized  by  a  geometric  imperfec¬ 
tion  of  the  form 


0  .0  ^  miTX  mry 

w  =  0  t  cos  — ^  cos  — -r- 

za  2b 


(15) 


where  6  is  a  measure  of  the  size  of  the  imperfection 
and  t  is  the  plate  thickness.  All  other  initial  stresses 
and  displacements  are  zero. 

=  u°=  0  (16) 

XX  yy  xy 

The  approximate  form  of  the  stresses  are  taken  to  be 
*  '  •  z  imrx  ntry 


•  •  z  miTX  n'lTv 

T  =  a,  —  cos  — r-  cos 
yy  2  t  2a  2b 


z  miTX  ,  niry 

T  =  a,  -  Sin  — sin 
xy  ^  t  2a  2b 


miTx 

2a 


(17) 


The  incremental  edge  displacements  can  be  specified  as 


u 

V 

w 


-  6*x 


=  0 


mTTx  nTry 

6t  cos  — cos 

2a  2b 


(18) 
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with  B  being  a  measure  of  the  size  of  the  edge  displacement. 
The  observer  may  notice  that  under  the  above  conditions,  the 
plate  is  allowed  to  move  in  such  a  way  that  midplane  stresses 
are  zero  in  shear  and  the  y-direction.  This  can  only  occur 
if  the  boundary  conditions  along  y  =  ±  b  are  held  by  a  very 
flexible  support.  In  order  to  create  v  =  0  displacement,  an 
external  source  is  necessary  to  relieve  any  restraining 
and  midplane  stresses.  With  this  artificial  boundary,  it 
is  possible  to  totally  isolate  the  plate's  imperfections  and 
remove  post  buckling  stresses  resulting  from  restrained 
boundaries . 

The  variational  principal 


61  =  0 


yields  for  this  problem 


310  gio 


which  is  a  system  of  differential  equations  in  the  unknown 

*  •  •  •  • 

incremental  parameters  6,  a,  Uj,  otj,  o;3. 

The  stress  and  displacement  distribution  can  be 
determined  by  adding  the  incremental  quantities  to  the 
previous  conditions. 


u  =  e  X 


v  =  0 

„  miTX  niry 

w  =  6t  cos  —jg  cos  — ^ 
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XX 


yy 


xy 


=  -a  +  Oj  ^  cos  miTx  cos  mry 


2a 

=  a,  z  cos  miTx  cos  miy 
t  2a  “2b 

=  Oj  £  sin  imrx  sin  niry 
t 


2b 


(21) 


2a  2b 

Consequently,  the  stress  and  displacement  at  any 
instant  of  time  can  be  determined  by  the  interative  solution 
of  a  set  of  five  algebraic  equations  in  6,  a,  ct^, 

Therefore,  the  problem  may  be  represented  in  matrix  format 


as 


a^j  (6,  a,  a^,  ,  a^) 


6 
a 
OL  . 


(22) 


with  the  initial  conditions  being 


6  =  6 


,0 

'  3 


Determination  of  a. .  Matrix 
- 13  - 

In  order  to  obtain  the  a^^  coefficients,  Eq ' s  (3,  17,  18) 
are  substituted  into  Eq.  (14)  with  the  required  partial  der¬ 
ivatives  of  Eq.  (18).  To  obtain  the  C^^'s,  Eq's  (10)  or  (11) 
are  substituted  into  Eq.  (14) .  Equation  (14)  can  be  developed 
for  the  elastic-plastic  case  by  substituting  Eq.  (21)  into 
Eq.  (5)  and  then  the  results  into  Eq.  (14) .  The  resulting 
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expression  is  then  introduced  into  Eq.  (13)  which  produces 
a  desired  form  of  the  functional  l°.  It  is  convenient  to 


take  the  variation  of  the  function  with  respect  to  the 
parameters  a,  a^,  ,  ot^  prior  to  the  required  integrations, 

as  several  terms  become  zero  in  the  derivative.  Switching 
the  order  of  variation  and  integration  is  possible  since  the 
integration  is  independent  of  the  variables  of  the  variation. 


A  detailed  derivation  and  listing  of  the  a^j  matrix  for 
elastic-plastic,  elastic  and  orthotropic  materials  is 
included  in  Appendix  B. 

At  this  point  the  problem  is  reduced  to  an  iterative 
solution  of  the  system  of  algebraic  equations  defined  by 
Eq.  (22) .  For  this  thesis  the  solution  to  the  matrix 
expression  uses  the  Fortran  IV  linear  equation  solver 
LEQTIF  from  the  IMSL  library.  The  solution  of  the  set  of 
linear  equation,  6,  o,  a^,  ,  a^,  is  added  to  the  previous 

values,  6,  o,  a,,  a  ,  a,  and  substituted  in  the  a . ^ 

12  3  A  J 
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equations,  resulting  in  new  values  for  aj^j  which  are  sub¬ 
sequently  used  in  the  next  iteration.  This  procedure  is 
continued  until  the  load-deformation  curve  is  determined. 

•k 

The  number  of  iterations  is  dependent  upon  the  size  of  B 
A  significant  advantage  of  this  procedure  is  that 
the  expressions  for  stress  and  displacement  allow  analyti 
integration  of  1°  which  eliminates  the  necessity  for  per¬ 
forming  expensive  numerical  integration  at  each  step  and 
permits  the  possibility  of  adapting  the  analysis  to  rel¬ 
atively  simple  computational  devices,  i.e.  advanced  hand 
calculators  or  desk  computers. 
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III.  Results  and  Discussion 


This  section  contains  a  suiranary  of  results  found  from 
the  different  applications  of  this  study.  A  computer  program 
was  created  which  utilized  the  analytical  technique  presented 
in  Section  II  to  calculate  the  load  versus  displacement  data 
for  the  various  plate  configurations  of  interest. 

Verification  of  Computer  Program 

Since  the  work  by  Neale  on  elastic-plastic  plates  (9) 
uses  the  same  analytical  technique  as  this  thesis,  his  results 
are  used  as  a  comparison,  to  verify  that  the  developed  computer 
program  produces  reliable  values. 

Sample  runs  were  made  with  square  plates  incorporating 
2a/t  ratios  of  20,  30  and  40.  The  results  are  compiled  in 
Table  I,  comparing  the  critical  compressive  loads.  The  critical 
compressive  load  is  taken  to  be  the  point  at  which  further 
increase  in  edge  displacement  results  in  a  decrease  in  the 
compressive  load,  a.  It  is  apparent  that  the  program  produces 
results  very  favorable  with  those  of  Neale.  In  fact,  all  the 
solutions  are  equal  to  or  lower  than  the  published  results.  As 
is  expected  with  a  Ritz  approximation,  the  solution  approaches 
convergence  from  above  and  consequently  the  results  that  are 
lower  than  Neale's  are  assumed  to  be  slightly  closer  to  exact 
values,  due  to  the  convergence  technique. 

Of  significant  interest  in  the  study  was  the  solution 
convergence  criteria.  To  arrive  at  a  solution. 
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a  particular  increment  size,  B  ,  was  chosen  and  the  program 
carried  out  the  iterations  until  the  stress  of  the  ith  iter¬ 
ation  was  less  than  the  stress  of  the  i-1  iteration.  This 
point  was  taken  as  the  relative  critical  stress.  The  step 
size  was  reduced  by  one-half  and  a  new  relative  critical 
stress  was  found.  The  process  continued  until  the  jth  rel¬ 
ative  critical  stress  fell  within  1%  of  the  j-1  relative 
critical  stress.  This  point  is  called  the  critical  com¬ 
pressive  stress  for  the  plate  (see  Figure  2) .  The  choice 

*  * 

of  B  is  of  course  important,  in  that  by  choosing  a  B  too 

large  one  obtains  a  curve  diverging  from  the  correct  solution. 
On  the  other  hand,  selecting  B*  too  small  required  many  more 
iterations  than  necessary  to  reach  convergence.  Since  one 
advantage  of  this  technique  is  the  simplicity  of  the  numer¬ 
ical  iteration  process,  the  requirement  for  more  iterations 
than  are  needed,  presents  a  drawback  to  adapting  the  procedure 
to  less  expensive  computing  devices. 

Figure  3  illustrates  a  typical  load  versus  displacement 
curve  for  an  elastic-plastic  plate  with  a  =  1.0,  b  =  1.0, 
t=  .1,  E=  30x10*,  v=  .3  and  Ramberg-Osgood  material 
constants,  from  Equation  (6),  of  E/E^  =  300  and  k  =  3.  The 
non-dimensional  load  cf/E^  is  plotted  against  the  displacement 
of  the  center  of  the  plate  in  the  z-direction.  Curves  are 
shown  for  initial  imperfections  of  1  half-sine  wave  in  the 
X  and  y  direction,  (i.e.,  m,  n  =  1) ,  of  magnitude  6°  =  10“^, 
10“ '*  and  10“ ®  inches. 


19 


j 


Stress 


Lateral  Displacement 

Fig.  2.  Convergence  to  Critical  Stress  through 
Repeated  Iteration 
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with  the  above  results  indicating  close  comparison 
with  published  data,  the  effect  of  varying  the  parameters 
was  evaluated. 

Elastic-Plastic  Materials 

Load  versus  displacement  curves  were  obtained  for 
plates  with  2a/t  ratios  of  10,  20,  40,  100,  200,  400.  The 
shape  of  the  individual  curves  are  consistent  with  those 
shown  in  Figure  3.  Increasing  initial  imperfection  size 
coincides  with  smaller  critical  stress  and  larger  displace¬ 
ment  at  collapse.  Since  the  shape  of  the  curves  do  not 
vary  greatly  from  those  shown,  individual  plots  are  not 
shown  for  each  case.  If  the  a  and  b  dimensions  of  the  plate 
are  increased,  keeping  a/b  equal  to  1,  the  stress  required  to 
reach  the  critical  stress  declines  rapidly  for  2a/t  less 
than  100.  The  rate  of  decrease  slows  and  approaches  zero 
as  2a/t  becomes  large,  (see  Figure  4) .  If  the  shape  of  the 
imperfection  is  changed  by  increasing  the  number  of  half¬ 
sine  waves,  the  critical  stress  is  increased  significantly. 
For  all  plate  sizes,  a/b  =  1,  the  minimum  critical  stress  is 
found  at  m,  n  =  1.  Further  discussion  of  the  effects  of  m 
and  n  is  found  later  in  the  text. 

If  the  ratio  of  length  to  width,  a/b,  is  varied,  a 
significant  difference  in  critical  stress  is  found  depending 
on  whether  the  length  is  held  constant,  (2a  =  constant) ,  and 
the  width  varied,  or  if  the  width  is  held  constant,  (2b  = 
constant) .  These  effects  are  shown  in  Figure  5  with  the 
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Plate  Size,  2a/t 

Fig.  4.  Effect  of  Plate  Size  on  Critical  Stress 
for  Elastic-Plastic  Plate 


aspect  ratio  plotted  against  the  normalized  critical  stress. 
It  should  be  noted  that  as  a/b  increases  the  critical  stress 
is  minimized  at  increasing  values  of  m.  For  a/b  less  than 
approximately  1  the  critical  stress  is  lower  when  the  width 
2b  is  constant,  whereas  for  a/b  greater  than  1  the  stress  is 
increasingly  lower  when  the  length,  2a,  is  constant.  This 
difference  can  be  explained  by  comparing  the  collapse  of  an 
elastic-plastic  plate  with  that  of  an  elastic  plate  and  will 
be  discussed  fully  in  the  following  section. 

Elastic  Material 

Figure  6  shows  the  stress  versus  displacement  curve 
for  an  elastic  plate  of  dimension  similar  to  those  of  the 
elastic-plastic  plate  in  Figure  3.  The  significant  diff¬ 
erence  lies  in  the  fact  that  the  three  curves  for  initial 
imperfections  of  6°  =  10~^,  10  ,  and  10”®  all  rise  asym- 

totically  to  the  same  value.  It  can  be  shown  that  this 
value  is  the  classical  bifurcation  stress  described  by  the 
following  equation  (Reference  3) ; 

0 =  kir^D  (24) 

(2b)A 

where  k  =/nvb  +  (25) 

\a  mb/ 

D  =  Et’  (26) 

12(l-v‘') 

A  =  area  of  load  application,  2bt 
Since  the  stress  functions  assumed  for  this  problem 
do  not  permit  evaluation  of  postbuckling  stresses,  the 
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•  Classical  Bifurcation 
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Fig.  6 


a  =  1.0 
b  =  1.0 
t  =  .1 


m,n  =  1 


- 1 - 1 - 

.1  .2 

Lateral  Displacement,  6 

Stress  vs.  Displacement  for  Typical 
Elastic  Material 


stress-displacement  curve  does  not  rise  above  the  bifurcation 

stress  as  would  be  found  in  an  actual  plate.  For  this  reason 

the  stress  associated  with  a  lateral  displacement  of  an  elastic 

plate  equal  to  the  thickness  of  the  plate  is  defined  as  the 

practical  critical  stress,  a  .  It  is  assumed  that  a  de- 

pr 

flection  oeyond  this  value  could  be  considered  as  beyond  the 
practical  plate  limits. 

The  shape  of  the  load-displacement  curves  for  the 
different  initial  imperfections  are  all  similar  to  those 
iii  Figure  6,  therefore  not  all  of  the  curves  are  shown. 

In  Figure  7,  load-displacement  curves  are  shown  for  two 
plates  of  dimension  2a  =  5-,  2b  =  5.  and  2a  =  10.0,  2b  =  10.0, 
considering  an  initial  imperfection  of  lO”"*.  The  elastic 
curves  are  compared  with  elastic-plastic  curves  for  the  same 
plates.  If,  in  Figure  7,  one  compares  the  critical  stresses 
of  tne  smaller  plate  with  those  of  the  larger  plate,  it  is 
apparent  that  the  effects  of  plasticity  are  significantly 
more  pronounced  in  the  smaller  plate  than  in  the  larger.  In 
fact,  the  difference  in  the  large  plate  is  only  approximately 
7%  as  compared  to  approximately  25%  for  the  smaller  plate. 
Comparing  the  whole  spectrum  of  2a/t  ratios.  Figure  8,  the 
effect  is  obvious.  It  is  clear  that  in  a  plate  requiring 
large  stress  to  reach  the  critical  stress,  the  effects  of 
plasticity  are  much  more  pronounced  and  that  in  a  plate 
requiring  small  stress  to  reach  critical  limits  the  plasticity 
effects  become  almost  negligible.  In  other  words,  the  thicker 
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Plate  Size,  2a/t 

Fig.  8,  Comparison  of  Elastic  and  Elastic-Plastic 

Critical  Stresses  as  a  Function  of  Plate  Size 
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the  plate  the  greater  the  plasticity  effects.  This  can  be 
seen  by  examining  the  stress-strain  curves  for  an  elastic 
and  elastic-plastic  material.  Figure  9. 

As  with  the  elastic-plastic  plates,  the  effect  of 
changing  the  a/b  ratio  was  computed  for  elastic  plates.  The 
results  are  shown  in  Figure  10  along  with  the  elastic-plastic 
results.  It  was  found  that  for  an  elastic  material,  it  did 
not  make  a  difference  which  dimension,  a  or  b,  was  kept 
constant  and  that  in  fact  the  same  curve  was  found  for  both 
cases . 

To  explain  the  fact  that  the  method  of  changing  the 
a/b  ratio  affected  the  outcome  in  the  elastic-plastic  case 
and  not  in  the  elastic  case,  one  must  look  at  the  effect  of 
the  m  and  n  values.  By  changing  the  shape  of  the  imperfection 
from  1  half-sine  to  3  half-sines  and,  likewise,  5  half-sines, 
the  effect  is  that  it  determines  the  collapse  mode,  that  is, 
the  final  shape  of  the  plate.  It  can  be  seen  from  Figure  10, 
that  as  the  a./ ^  ratio  increases  the  minimum  stress  for  collapse 
occurs  at  higher  m  values  and  that  for  an  elastic  material  the 
mimimum  values  become  asymtotic  to  some  mimimum  value.  From 
this,  it  can  be  seen  that  a  plate  with  an  initial  imperfection 
of  1  half-sine  wave  will  collapse  at  the  same  stress  as  a  plate 
whose  dimensions  are  odd  integers  larger  than  the  single  half¬ 
sine  plate,  provided  an  imperfection  exists  of  the  same  odd 
integer  shape.  For  example,  in  the  case  of  a  square  plate 
15  inches  by  15  inches  with  an  initial  imperfection  of  3 
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Stress-Strain  Relation  for  Elastic  and 
Elastic-Plastic  Material 
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Aspect  Ratio  a/b 
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half-sine  waves  in  the  x  and  y  directions,  the  collapse 
stress  is  the  same  as  a  plate  5  inches  by  5  inches,  with 
1  half-sine  imperfection  in  the  x  and  y  direction.  Con¬ 
sequently  a  plate  with  1  half-sine  wave  in  each  direction 
can  be  thought  of  as  an  "elemental  component"  for  larger 
plates.  (See  Figure  11). 

With  this  concept  in  mind,  one  can  justify  the 
difference  in  the  two  methods  of  changing  the  a/b  ratio  for 
an  elastic-plastic  plate.  Figure  12  illustrates  the  variation 
in  a/b  using  both  the  constant  a  and  constant  b  methods.  In 
the  b  =  constant  case,  for  a/b  =  3,  Figure  12f,  the  plate  is 
three  times  as  long  as  the  a/b  =  1  case,  Figure  12e,  the 
minimum  collapse  stress  is  found  when  an  initial  imperfection 
of  m  =  3  is  used.  Consequently,  the  "elemental  component"  is 
the  same  size  as  the  a/b  =  1  plate  and,  therefore,  the  critical 
stress  is  the  same  as  the  a/b  =  1  case.  Similar  conditions 
apply  for  m  =  5.  For  the  other  method,  (a  =  constant),  when 
a/b  =  3,  Figure  12c,  the  size  of  the  "elemental  component"  is 
1/3  as  large  as  the  a/b  =  1  condition.  Figure  12b,  and  as  was 
shown  in  Figure  8,  as  2  a/t  decreases,  the  plasticity  effects 
become  more  pronounced.  Therefore,  the  curve  for  the  elastic- 
plastic  with  a  =  constant.  Figure  10,  is  expected  to  decrease 
as  a/b  increases.  In  effect,  the  results  indicate  that  a 
plate,  whether  elastic  or  elastic-plastic,  with  initial 
imperfections  in  the  x  and/or  y  direction  will  collapse  at 
the  same  stress  as  the  "elemental  component",  that  is,  the 
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a  =  constant 


a/b  <  1 


a/b  =1  a/b  >  1 
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b  =  constant 


a/b  <  1 
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d. 
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Fig.  12.  Methods  of  Changing  a/b  Ratio 


integral  part  of  the  plate  with  a  single  half-sine  wave  in 
the  X  and  y  direction. 


Orthotropic  Materials 

The  orthotropic  plate  considered  in  this  thesis 
possesses  elastic  properties  along  the  major  material  axis. 
Consequently,  the  results  found  are  similar  to  those  of  the 
previous  section.  A  significant  interest  was  in  the  effect 
that  the  orientation  of  the  material  had  on  the  critical 
stress.  The  material  properties  are  considered  to  be 
similar  to  those  of  boron-epoxy  composite  with  Ej/E^  =  10, 
Gjj/Ej  =  .03,  Vj2  =  .3. 

It  was  determined  early  in  the  study  that  results  for 
off  axis  orientation  were  significantly  lower  than  expected. 
The  reason  for  this  can  be  seen  if  one  examines  the  diff¬ 
erential  equation  governing  an  orthotropic  plate  of 
arbitrary  orientation  (11) 
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where  the  subscripts  denote  differentiation  with  respect  to 
the  X  and  y  directions.  The  bending  stiffnesses,  ,  are 
described  in  detail  in  Jones  (11) .  The  major  difference 
between  this  differential  equation  and  that  for  an  ortho¬ 
tropic  plate  oriented  at  0°  and  90°  to  the  major  material 
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axes,  referred  to  as  specially  orthotropic,  is  in  the 

inclusion  of  the  twist  coupling  stiffnesses,  D  and  D 

16  2  6 

It  can  be  shown  that  D  and  D, ,  are  dependent  upon  the 
constitutive  coefficients  Cjj  and  3  of  Eq.  f4). 

However,  the  form  of  the  stress  function  assumed,  is 
such  that  these  cross  coupling  terms  do  not  carry  through 
the  integration  and  variation  of  the  functional,  illustrated 
in  Appendix  B.  Consequently,  the  stresses  at  an  orientation 
other  than  0°  and  90°  are  reduced  by  an  amount  equivalent 
to  the  cross  coupling  terms.  For  an  orientation  of  0°  and 
90°  these  terms  and  Dje'  found  in  the  differ¬ 

ential  equation  and  consequently  reliable  results  are 
expected . 

Mandell  (15)  made  an  experimental  investigation  on 
the  buckling  of  anisotropic  fiber  reinforced  plates  with 
boundary  conditions  similar  to  those  imposed  in  this 
report.  A  boron-epoxy  plate  11  inches  by  9.75  inches 
by  .096  inches  thick,  simply-supported  along  all  edges 
with  an  in-plane  displacement  allowed  in  the  normal  direc¬ 
tion  and  loaded  at  0  =  0°,  was  considered.  Results  are 
shown  in  Figure  13.  Comparison  is  made  with  curves  obtained 
using  an  initial  imperfection  of  10~^  and  10“^.  Note  that 
the  two  curves  provide  upper  and  lower  bounds  to  the 
experimental  results.  This  indicates  that  the  inclusion 
of  geometric  imperfections  does  indeed  result  in  load- 
deflection  curves  which  are  representative  of  actual  plates. 
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Fig.  13.  Comparison  of  Experimental  with 

Analytical  Results  for  Orthotropic 
Plate 


at  least  in  the  initial  stages  of  deflection.  The  critical 
stresses  are  shown  for  6=0°  and  90°  along  with  experimental 
values  for  other  angles  in  Figure  14.  A  Rayleigh-Ritz  theo¬ 
retical  solution  is  also  shown.  The  thesis  results  compare 
favorably  at  6  =  0° .  At  9  =  90°  the  results  are  considerably 
above  experimental  values  and  even  more  above  the  Rayleigh- 
Ritz  solution.  The  classical  buckling  mode  for  this  condition 
is  with  m  =  2;  however,  the  initial  imperfection  shape,  in 
this  case  m  =  3,  establishes  the  mode  at  collapse  and  con¬ 
sequently  the  critical  stress  obtained  at  90°  is  higher  than 
the  classical  value. 

For  an  a/b  ratio  of  1,  the  affect  of  2a/t  is  shown 
on  Figure  15.  For  this  case  both  the  9=0°  and  the  6  =  90° 
curves  are  coincident.  It  can  be  seen  in  Figure  16  that  this 
is  not  necessarily  true  for  other  values  of  a/b.  In  fact,  for 
a/b  less  than  1  the  critical  stress  for  6  =  90°  becomes 
significantly  less  than  the  6=0°  orientation.  In  Figure  16 
the  Ej/E^  ratio  is  10,  indicating  that  the  strength  in  the 
2-direction  is  much  less  than  in  the  1-direction.  If  one 
decreases  this  ratio,  as  with  a  boron/aluminum  composite, 

Ej/Ej  =  1.44,  the  difference  between  6=0°  and  90°  becomes 
less  pronounced.  (See  Figure  17.)  This  is  as  expected  since 
as  Ej/E^  tends  to  1,  essentially  isotropic,  the  curves  would 
merge  into  one.  The  results  found  for  the  boron/aluminum 
composite  compare  very  well  with  results  published  by 
Chamis  (14)  for  a/b  less  than  1.  For  a/b  values  greater 
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Fig. 14  Comparison  of  Orthotropic  Results 


Stress 


than  1, the  effect  of  the  imperfection  shape,  m,  is  more 
pronounced  but  approaches  the  same  minimum  stress  as 
published  by  Chamis  as  a/b  increases. 


IV.  Conclusions 


A.  Elastic-Plastic 

1.  The  size  of  an  imperfection  greatly  affects  the 
load  displacement  curve.  For  an  elastic-plastic  material, 
not  only  the  degree  of  displacement  at  any  instant,  but  the 
stress  at  which  collapse  occurs  is  dependent  upon  the 
imperfection.  A  plate  with  a  large  initial  imperfection 
will  collapse  at  a  lower  stress  and  larger  displacement 
than  a  similar  plate  with  small  initial  imperfections. 

2.  As  the  aspect  ratio  increases,  the  normalized 
stress  for  a  particular  m  value  decreases  to  a  minimum 
and  then  begins  to  increase  again.  The  minimum  stress 
required  for  collapse  is  found  at  increasing  values  of 

m  for  increasing  a/b  ratios.  In  all  cases  the  minimum 
stress  is  fou.. .  with  n  =  1. 

3.  The  shape  of  the  initial  imperfections  determined 
the  shape  of  the  plate  at  collapse.  That  is,  if  a  plate 
has  an  initial  imperfection  of  m  half-sines  in  the  x 
direction  and  n  half-sines  in  the  y  direction,  it  will 
collapse  in  that  mode,  regardless  of  plate  size  or  aspect 
ratio. 

4.  Plasticity  affects  become  more  and  more  signif¬ 
icant  as  the  effective  thickness  of  the  plate  increases. 

5.  A  plate  with  an  initial  imperfection  of  m,  n  =  1 
is  the  elemental  component  from  which  the  critical  stress 
of  plates  with  other  imperfection  modes  may  be  determined. 


B.  Elastic 


1.  The  affect  of  imperfection  size  is  similar  to  the 
elastic-plastic  material  except  that  the  collapse  stress 
becomes  asymtotic  to  the  classical  bifurcation  stress. 

Because  of  the  form  of  the  assumed  stress  functions,  no 
postbuckling  stresses  can  be  obtained.  Consequently,  the 
stress  at  a  displacement  equal  to  the  thickness  of  the 
plate  was  assumed  to  be  the  practical  critical  stress. 

2.  The  effect  of  the  magnitude  of  the  imperfection 
on  the  practical  critical  stress  becomes  more  pronounced 
as  the  imperfection  size  increases,  that  is, the  practical 
critical  stress  decreases  with  increasing  imperfection  size. 

3.  Critical  values  computed  are  dependent  upon  the 
imperfection  shape,  m.  Individual  critical  stresses  may  be 
higher  than  classical  values  because  not  all  mode  shapes 
are  considered. 

4.  The  effects  of  plate  size,  aspect  ratio  and 
imperfection  shape  are  all  similar  to  those  of  the  elastic- 
plastic  case. 

C,  Orthotropic 

1.  Imperfection  size  and  shape,  plate  size  and  aspect 
ratio  all  have  affects  similar  to  those  of  an  elastic  iso¬ 
tropic  material. 

2.  The  orientation  of  the  plate  has  a  significant 
effect  on  the  collapse  stress  at  a/b  ratios  less  than  1. 

The  relative  difference  between  Ej  and  E^  determines  the 
magnitude  of  the  difference. 
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Appendix  A 

Material  Properties 

1)  Elastic  Plastic 

E  =  30  X  10®  Ib/in^ 

Eq  =  100,000  Ib/in^ 

G  =  27/4E^ 

0 

V  =  0.3 
k  =  3 

2)  Elastic 

E  =  30  X  10®  Ib/in^ 

V  =  0.3 

3)  Orthotropic 

Boron-Epoxy 

Ej  =  30  X  10®  Ib/in^ 

E^  =  3.0  X  10®  Ib/in^ 

=  1-0  X  10®  Ib/in^ 

^2  =  0-3 

Boron-Aluminum 

Ej  =  35  X  10®  Ib/in^ 

E^  =  24.3  X  10®  Ib/in^ 
Gj2  =  11-6  X  10®  Ib/in^ 
=  0.24 


( 
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Appendix  B 

Calculation  of  a. .  Matrix 
ID 


Beginning  with  the  functional, 
t/2  b  a 

=  /  /  / 

-t/2  -b  -i 


where  t..E..  =  t  E  +2t  E  +  E 

ID  ID  XX  XX  xy  xy  yy  yy 


dxdydz 

(Bl) 

(B2) 


/3w 

_  9^ 

9w 

/  aw 

(ax  y 

ay  ^yy 

^ij  '^k,i  '^k,j  ^xx 

+  2C,,t^,t  +  2C,  T  ^ 

2  3  xy  yy  3  3  xy 

To  get  Iq  in  terms  of  known  quantities,  substitute  the 
assumed  functions 


(B4) 


•  ★ 


u  = 

-6  X 

• 

(B5a) 

« 

V  = 

0 

(B5b) 

« 

w  = 

« 

6t  cos 

mux 

“2i 

nuy 

2b 

(B5c) 

• 

T 

XX 

• 

=  -a 

+ 

•  z 

cos 

mux 

~~7a 

cos 

nuy 

2b 

(B6a) 

>1 

• 

=  “2 

z 

t 

mux 
cos  -Jg 

cos 

2b 

(B6b) 

• 

T 

xy 

=  i, 

z 

t 

mux 
sin  2a 

sin 

2b 

(B6c) 

u  = 

-6  X 

(B7a) 

V  = 

0 

{B7b) 

w  = 

6t  cos 

mux 

2a 

niTy 

2b 

(B7c) 

^xx 

=  -0 

+ 

z  mux 

“1  Z  -JE 

n^y 

2b 

(B8a) 
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(B8b) 


T  =  a„  z  cos  mirx  cos  mT\ 

yy  ^  t 


itnrx  cos  mry 
2a  “2E 


r  =  a.  z  sin  imrx  sin  mry 
^  t  -JE  -2E 


The  strain  rate  equations  are; 


XX 


•  "k 


3x 


9x 


9u 

-  z 

92w 

+  9w 

9w 

9x 

3x^ 

¥x 

9x 

9v 

z 

9  2w 

+  9w 

9w 

9y 

9y^ 

Ey 

9y 

%  . 

9u  + 

9v 

-  2z 

+  9w 

9w 

9y 

9x 

9x9y 

9y 

+ 

9w 

9w 

9x 

9y 

be 

found  from  Equations 

(B5 

-  B7) 

'  0 

9w  = 

-miT 

♦ 

6t  sin  mux  sin  nuy 

9x 

2a 

2a 

2b 

1, 

9w 

=  — : 

mr  St 

cos 

mux 

sin 

mry 

9y 

2b 

2a 

2b 

9y  9y 

9 =  0  6t  cos  tmrx  cos  niry 

3X^  4a^  2a  2b 

9^w  =  -n^TT^  6t  cos  mtrx  cos  niry 
45^  2a  2b 

9  w  =  mmr'^  6t  sin  imrx  sin  mry 


9x3 


xoy 


4ab 


2a 


2b 


9w  =  -miT  6t  sin  ttnrx  cos  mry 
9x  2a  2b 

9w  =  -rnr  8t  cos  imrx  sin  niry 
9y  2E  ~2a  2b 

If  Equations  (BIO  -  16)  are  substituted  into  {B9) ,  the 
results  become 


•  * 


=  -6  -  z  -m^TT*  6t  cos  mux  cos  mry 

A  -  V  — Trr-  — 


XX 


“2a 


2b 


(B8c) 

(B9a) 

(B9b) 

{B9c) 

(BIO) 

(BID 

{B12) 

(B13) 

(B14) 

(B15) 

(B16) 


■ 


xy 


-mir 

6 

t 

sin 

mirx 

cos 

mry  ' 

\/-mir  6 

t 

sin 

mirx 

cos 

mry 

\ 

2a 

2a 

2b  , 

A  2a 

2a 

2b 

= 

- 

■z 

( 

-n 

• 

6  t  COS  mirx  cos  mry 

1 

\ 

2a 

2b 

) 

-mr 

5 

t 

COS 

mirx 

sin 

niz 

V-niT  6 

t 

cos 

mirx 

sin 

niz 

\ 

2b 

2a 

2b 

A  2b 

2a 

2b 

= 

f-2z 

/  mmr^ 

6  t 

sin 

mirx 

sin  mry  '\ 

L 

4ab 

2a 

2b  J 

+  ( 

^-mir 

6 

t 

sin 

mirx 

cos 

mr^' 

\/-niT  6 

t 

cos 

mirx 

sin 

mry 

b  2a 

2a 

2b, 

A  2b 

2a 

2b 

+  (' 

’  -mir 

• 

6 

t 

sin 

mirx 

cos 

mry  ^ 

kZ-mr  6 

t 

cos 

mirx 

sin 

mry 

V 

2a 

2a 

2b, 

A  2b 

2a 

2b 

)] 


Equation  (B3)  can  now  be  rewritten  with  the  use  of  equations 
(B8)  and  (BIO  -  16)  . 


ID 


-i  V.  .  2  cos  imrx  cos  nTry  \ 

[\  ‘  t  ~2a  2b  J 

•  • 

*f  -itnr  5  t  sin  mvx  cos  mry  \/-ni'JT  6  t  sin  mirx  cos  niry  \ 

\  2a  2a  2b  A  2a  2a  2b  J 

+  2/a,  z  sin  mivx  sin  mryV-imr  6  t  sin  imrx  cos  mry  \ 

V  '  t  Ti  ~2b  A  2i:  -JE  ) 

•  /-mr  5  t  cos  imrx  sin  mry  \ 

V  2b  2a  2b  / 

• 

+/a  z  cos  mrrx  cos  nny  \/~n  6  t  cos  irnrx  sin  mry  \ 

\  t  2a  2b  A  2b  2a  2b  / 

•  , 

•/-iw  5  t  cos  mirx  sin  mry  \ 

\  2b  2a  2b  J 

Similarly,  the  substitution  of  Equation  (B6,  8)  into 

Equation  (B4)  results  in: 

W(t)  *^|c.,/-a  +  a,  z  cos  mirx  cos  mry  \ 

I  ‘‘V  '  t  -JE  ~iB) 

,/-o  +  a  z  cos  mirx  cos  mry  \ 

V  ‘  t  "IE  ) 


(B17a) 


(B17b) 


{B17c) 


(B18) 
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•  • 

-a  +a  z  cos  imrx  cos  niry  Vg.  £  cos  m-rrx  cos  mry  \ 
*  t  2a  2b  A  t  2a  2b  / 


•  • 

+  C  /a,  2  cos  mTTx  cos  nTry  \/ot,  z  cos  ihttx  cos  mry  \ 

“I  t  2a  2b  A  t  2a  2b  j 

+  2C,  -a  +  a,  2  cos  mTrx  cos  n7ry\/ot,  z  sin  mTrx  sin  nTry  \ 

'  t  “2i  ^A  '  t  -2i  it; 

+  2C.  ,/a  2  sin  m7fx  sin  mry\/a,  z  cos  imrx  cos  niry 'N 

(  ^  t  2a  1b )(  t  2a  it  j 

+  2C,,/a,  2  sin  irnrx  sin  mry  \/a,  z  sin  mirx  sin  niry  N]  (Bl9) 

\  t  2a  it  ;V  t  2a  it  j  ] 

The  expansion  of  Equation  (B2) ,  using  (B6)  and  (B17) ,  produces 


C 


T  E  =a6*  +  a2  /  6  t  cos  imrx  cos  niry  \ 

XX  XX  ^  4a'^  2a  2b  ; 

-/a  6  6  sin^  imrx  cos^  mry  \ 

V  ~4a^  2a  2b  / 

“6*a  2  cos  mrrx  cos  niry  +  a  z^m^TT^  6  cos^  irnrx  cos*  nrry 

t  2a  2b  4a2  2a  2b 

+a,  z  in*7r*t  6  5/cos  mrrx  cos  nrryVsin*  mrrx  cos*  mTy\(B20a) 

'  4^  \  -JK  it;i  n  it; 

•  •  •  • 

2t  E  =  -2a,  z*  6  sin*  mrrx  sin*  mry 

xy  xy  '  H  lt 

+2  a  zt  mniT*  6  6  sin*  mrrx  sin*  mry  cos  mrrx  cos  mry  (B20b) 

4ab  2a  2b  2a  2b 

T  E  =  a,  z*n*rr*  6  cos*  mrrx  cos*  nrry 

yy  yy  Tb~  li  it 


+  02 


zn  *  rr  *  t  <5 
4b  2 


6 


cos  * 


mrrx  cos  nrry  sin*  nrry 
2a  2b  2b 


{B20c) 


Since  not  all  the  terms  of  the  integrand  of  1°  involve  all 


the  incremental  parcimeters,  it  simplifies  the  mathematics  to 

•  •  • 

take  the  variation  of  the  integrand  with  respect  to  6 ,  a ,  a,, 
a^,  Oj  before  performing  the  integration,  with  the  result  that 


( 


c 


several  of  the  partial  derivatives  become  zero.  Switching 
the  order  of  the  integration  and  the  variation  does  not 
affect  the  solution  since  the  variation  is  taken  with 
respect  to  parameters  independent  of  the  integration 
parameters  x,  y,  and  z.  Consequently, 


(B21) 


with  similar  expressions  for  the  remaining  incremental  param¬ 
eters,  Carrying  out  the  above  actions  involves  considerable 
algebra  and  simple  calculus  and  therefore  all  the  required 
steps  are  not  shown.  It  should  be  observed  that  by  assuming 
a  function  for  stress  intensity  of  the  same  form  as  the 
imperfection,  one  obtains  relations  which  reduce  the  com¬ 
plexity  of  the  mathematics.  This  however,  eliminates  the 
possibility  of  relating  the  stress  boundary  to  a  post 
buckling  consideration.  Thus,  one  of  the  drawbacks  to  a 
stress  intensity  function  along  a  boundary  is  to  produce 
what  may  be  thought  of  as  unrealistic  displacement  conditions. 
The  approach  followed  herein  is  an  attempt  to  isolate 
imperfections  and  as  such  the  stress  expressions  assumed, 
can  be  justified. 
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( 


c 


5 

! 

< 

1 


The  W(t)  term  of  the  functional  1°,  contains  the 
constitutive  coefficients  C..,  that  vary  with  material 
property.  To  illustrate,  the  calculation  of  a^ ^  is 
shown  as  an  example.  From  Equation  (23)  of  Section  II 
it  can  be  seen  that  a  is  the  coefficient  of  ct,  in  the 
variation  of  1°  with  respect  to  a.  Upon  taking  the 
variation  of  the  integrand  of  1°  with  respect  to  a,  the 
expressions  become; 


9(t. .  E. .) 

_ ^3  ^2- 

3a 


•  4r 


e 


z  m^ir^ 


4a‘ 


6t  cos  miTx  cos  niry 
2a  2b 


-m^Tf^t^  6  6  sin^  mirx  cos^  mry 
4a^  2a  2b 


(B22) 


3  ( T  .  .  V,  .V,  . )  =  0 

13  k,i  k,]' 

3a 


3  W(t)  =  -C  a  z  cos  miTx  cos  nir 
3*  '  t  2a  2 


(B23a) 


The  only  term  containing  a  coefficient  of  is  Equation 
(23a),  consequently  to  find  ^  one  must  integrate  (B23a) . 
However,  must  be  expressed  in  terms  of  the  material 

constitutive  relations.  For  an  elastic-plastic  Equation 
(B23a)  is  expanded  to; 


3  (W(t)  =  - 
3a 
,2 


1  a.  z  cofa  mTTX  cos  ivir 
E  t  “2a  2 


+  G  /  4  a  cos  mirx  cos  niry  -  4  a  a  a  z 
9  V  2a  2b  ^  t 


cos^  niry 
2b 


cos'^  mirx 
2a 
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+2  a  ct,  a,  cos  irnrx  cos  mry  -  4  a  a,  a,  z^  cos^ 

"  '  -2E  '  ' 

imrx  cos^  niry 
2a  2b 

+a  4  z^  cos^  itnrx  cos^  mry  -  2a, a, a,  z^  cos  ^  mrrx 

2a  2b  2a 

cos  ^  mry 
2b 

+2a  a,  a  z^  cos^  mrrx  cos^  mry  -  2a  a  a  z^  cos^ 

2a  2b  ^ 

mrrx  cos^  mry 
2a  2b 

+af  a,  2^  cos  ^  mrrx  cos  ^  mry  \] 

The  above  expression  is  integrated  as  shown 


sf  W  (t  )  ^  dxdydz 

90 


During  the  integration  with  respect  to  z ,  all  terms  con¬ 
taining  odd  powers  of  z  go  to  zero  due  to  the  symmetry  of 
the  plate.  Consequently,  the  expression  is  reduced  to: 
t/2  b  a 


/  /[« 


'8aa .  a , 


4aa  01,^2^  cos^  mrrx 
^  Vt  2a 


•'t/2  -b  -a 


•  cos^ 


nrry 

■2b 


dxdydz 


Factoring  constants  outside  the  integral  and  integrating 
with  respect  to  z  produces: 


(B23b) 


{B23c) 
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Gt 

TM 


aajO  j 


imrx  cos 
2a 


mry  \  dxdy 
2b  )  (23d) 


Integrating  in  x  and  y,  using 


/cos^  CA  dX  =  _X  +  sin  2CX 

2  4C 


where  mj:  =  C 
2a 

For  odd  values  of  m  and  n,  the  expression  reduces  to: 


81°  a.  =  Gc^2aj-a2  ^ctj  (tba)  (23e) 

8d  '  27  V  I 

Since  the  only  values  of  m  and  n  that  will  fulfill  the 
boundary  conditions  are  odd  integers,  i.e.,  1,  3,  5..., 
Equation  {B23e)  is  the  complete  form  of  the  expression. 

Upon  completion  of  all  calculations  involving  the  coeffi¬ 
cient  of  the  a2j  row,  the  quantity  in  parentheses,  (tba) , 
may  be  divided  out  of  the  equation,  leaving: 


2  3 


1_  Ga(2a,-a^) 


(B24) 


In  a  similar  manner  the  remaining  a.  terms  can  be 
found.  For  an  elastic-plastic  material,  the  a.,  coefficients 
are  calculated  to  be: 


12 
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a  =  2  Gaa3 

52  ^ 

a  5  3  =  G  01 

320 

a  =  G  a  f  2a.  -  a  .\, 

“  320  '  ^  / 

a.,  =  1+v  +  27  Ga^  (B25) 

”  -X"  Tm  ^ 

If  one  applies  Equation  (B23a)  to  an  elastic  material 
and  makes  the  substitution  for  Cjj  the  calculation  simplifies 
as  follows 


3  W(t) 

90' 

where 


z  cos  mTTX  cos  n7ry 
t  2a  2b 


{B23a) 


Again,  the  integral  reduces  to  zero  when  integrated  with 
respect  to  z.  Consequently, 


a^3  =  0 


Likewise,  the  remaining  a^^  coefficients  for  a  material 
displaying  elastic  properties  are  determined  as: 


a 


1 1 


o , 


59 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


1  5 

2  2 

2  4 

3  1 

3  3 

3  5 

4  2 

4  4 

5  1 

5  3 

5  5 


^2  1 

/mirt 

6  V  "'t)  /  » 

Z  1 

= 

4  c,, 

^2  3 

= 

0 

0 

0 

2  5 

/mirt 

^32 

= 

0 

\  2a  J 

= 

1 

o 

^34 

= 

0 

^4  1 

/nirt 

(  2b} 

= 

0 

a,  , 

a 

4  3 

3  4 

* 

a 

, 

0 

2  2 

4  S 

/mirt  X/nTTt  \ 

0 

V,  2aJ\  2h) 

5  2 

0 

^5  4 

= 

0 

(B26) 


where  is  given  by  Equation  (10)  of  Section  II. 

Since  the  only  difference  between  the  elastic  and 
orthotropic  constitutive  relations  are  through  the  inclusion 
of  the  differences  in  properties  in  the  longitudinal  and 
transverse  direction,  and  in  the  coupling  of  the  properties, 
very  little  difference  is  found  in  the  calculation  of  the 
a. .  matrix.  Consequently,  the  a. ,  terms  for  a  material 
displaying  orthotropic  properties  can  be  expressed  as  those 
found  in  (B26)  where  the  C..'s  are  given  by  Equation  (11) 


of  Section  II. 
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